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Heat  Capacity  of  Water  at  Extremes  of  Supercooling  and  Superheating 


C.  A.  Angell,  W.  J.  Sichina  and  M.  Oguni 
Department  of  Chemistry 
Purdue  University 
West  Lafayette,  Indiana  47907 


Abstract 


;A 


\/>Low  temperature  heat  capacities  for  water  and  lbjCi  have  been  remeasured 

■'V 

using  different  sample  preparation  and  data  analysis  procedures  from  those 
used  in  the  original  studies  of  this  anomalous  region.  Self-consistent  re¬ 
sults  have  been  obtained  with  reproducibilities  of  1%.  Comparisons  are  made 
with  the  behavior  in  the  superheated  region  at  1  atm  pressure  calculated 
from  a  refined  equation  of  state  by  Haar  et  al.  It  appears  that  the  kinetic 
metastability  limit  of  the  liquid  (nucleation  terminated)  at  both  hot  and  cold 


extremes  may  be  reached  at  the  same  value  of  the  liquid  heat  capacity 


Y 
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1°  1972  Rasmussen  et  al^^  reported  that  the  heat  capacity  of  water  increases 
very  rapidly  on  supercooling  below  ^  -15°C  and  reaches  values  of  1.5  cal  gm-1deg 
(113  J  mol  ^deg  at  -38°C  at  which  temperature  crystallization  of  ice  I  termi¬ 
nated  the  measurements.  At  the  other  extreme,  equally  striking  increases  in 
heat  capacity  of  water  near  the  superheating  limit  have  been  determined  by 

extrapolations  to  low  pressure  (based  on  corresponding  states  arguments)  of  high 

(2  3) 

pressure  stable  range  data.  ’  It  had  earlier  been  suggested  by  Anasimov 

(4) 

et  al.  that  water  would  show  causally  related  anomalies  in  Cy  at  high  and 
low  temperature  extremes. 

The  purpose  of  this  report  is  to  present  new  and  improved  data  for  the 
low  temperature  range  and  to  point  out  that  the  phase  change  which  terminates 
the  metastable  region  at  each  extreme  occurs  in  the  region  of  sharply  increasing 
heat  capacity. 

There  is  a  particular  need  for  improved  data  at  low  temperatures.  The 
results  of  Rasmussen  et  al,^  which  have  helped  generate  much  of  the  current 
interest  in  supercooled  water,  were  essentially  of  exploratory  character.  They 
revealed  an  interesting  phenomenon  but  exhibited  differences  of  up  to  10% 
between  data  sets  from  different  methods  of  measurement,  and  different  instru¬ 
ments.  Nevertheless,  for  want  of  better,  these  data  have  been  used  as  a  basis 
for  predicting  and  analyzing  the  behavior  of  other  properties,  such  as  sound 
velocities, possibly  to  the  detriment  of  comparisons  with  experiment. 

For  instance,  if  the  values  of  Cp  used  were  to  contain  an  error  of  +  5%  growing 

in  between  -25  and  -33°C  in  the  region  of  fast-increasing  C  ,  then  a  predicted 

P 


2 


velocity  minimum  based  on 


would  disappear,  if  were  assumed  constant.  If  were  (more  correctly) 

subjected  to  the  same  correction,  then  the  predicted  velocity  minimum,  and  the 

associated  Kg  maximum,  would  be  displaced  to  lower  temperature,  and  a  phenomeno- 

(7  8) 

logical  relationship  to  the  maximum  in  Cv,  *  which  would  now  be  revealed 
more  clearly  at  'v  -20°C,  would  become  more  obvious.  Likewise,  the  thermodynamic 
absurdity 

0  >  Cp/Cv  -  [1  -  Ta2/pCp<T]_1  (2) 

obtained^  by  boldly  extrapolating,  to  T  <  -40°C,  power  law  equations  fitting 
a,  Cp,  and  tcT  data  measured  at  higher  temperatures,  will  not  appear. 

We  have  therefore  carefully  repeated  the  low  temperature  measurements  a 
number  of  times  using  differential  scanning  calorimetry  of  dispersed  droplet 
samples,  prepared  and  measured  by  different  investigators,  using  different  sample 
preparation  techniques  and  somewhat  different  data  analysis  methods.  We  pre¬ 
sent  here  sets  of  emulsion  data  which  are  self-consistent  to  ±  1%,  and  agree 
to  0.5%  with  bulk  sample  data  to  the  low  temperature  limit  of  the  latter,  -13°C. 
At  the  low  temperature  extreme,  the  new  data  differ  by  almost  5%  from  the 
average  values  previously  reported. 

EXPERIMENTAL 

In  the  low  temperature  emulsion  measurements,  the  water  contribution  to 
the  total  heat  capacity  of  the  two  phase  sample  has  been  separated  not  only  by 
the  technique  of  comparing  frozen  and  unfrozen  sample  data  used  previously^ 


but  also  by  the  technique  of  off-setting  the  non-aqueous  component  with  an 
equal  mass  of  dispersant-surfactant  solution  in  the  reference  pan.  Comparison 
of  results  with  literature  data  in  the  temperature  range  above  0°C  suggest 
that  these  techniques  are  of  approximately  the  same  overall  reliability.  The 
former  suffers  from  the  addition  of  the  random  errors  from  two  scans,  while  the 
latter  suffers  from  a  combination  of  errors  from  two  weighings  of  small  quanti¬ 
ties,  and  additional  uncertainties  from  possible  changes  in  emulsion  volume 
fraction  of  water  during  transfer  to  the  sample  pan. 

Measurements  were  performed  using  the  Perkin  Elmer  DSC-.2  instrument 
equipped  with  the  P.E.  Scanning  Auto  Zero  Accessory.  Samples  of  mass  ~  15  mgm 
were  down-scanned  from  *v  20°C  at  10  deg  min  *  on  a  sensitivity  setting  of 
2  >  mcal/s.  In  the  case  of  investigator  I's  determinations,  the  pen  displace¬ 
ment  with  the  instrument  in  the  isothermal  hold  condition  was  carefully  checked, 
for  both  empty  and  loaded  pans,  at  short  intervals  between  the  initial  and 
final  temperatures  of  the  scanning  run.  Under  favorable  instrument  conditions 
this  is  unnecessary.  However,  in  practice,  variations  in  these  relative  iso¬ 
thermal  hold  pen  displacements  which  would  imply  heat  capacity  errors  of  up 
to  5 X  for  small  samples  studied  in  uncalibrated  wide  temperature  range  scans, 

can  at  times  be  observed,  and  this  practice  is  necessary  to  obtain  the  best 

(q) 

performance  from  the  instrument. 

The  question  of  whether  the  lowest  temperature  results  are  influenced 
by  partial  crystallization  during  cooling  cannot  be  satisfactorily  answered 
by  the  obvious  procedure  of  making  measurements  during  temperature  increase 
rather  than  decrease,  though  data  above  -33°C  were  confirmed  by  this  method 

(lc) 

in  earlier  work.  The  reason  is  that  any  measuements  must  start  from  an 

"isothermal  hold"  position,  and  a  modest  time,  hence  temperature,  interval 
( 'y  2°C)  must  be  allowed  for  the  instrument-sample  ensemble  to  reach  its  correct 
steady-state  condition  for  measurement.  The  isothermal  hold  temperature  cannot  be 


4 


set  too  close  to  T„  because  the  total  elapsed  time  near  T„  determines  the 
n  H 

probability  of  nucleation.  Thus  the  temperature  is  already  above  the  region, 

T  <  35°C,  in  which  the  threat  of  incipient  crystallization  influence  is  real, 
before  the  measurement  can  commence.  A  better  confirmation  of  the  absence 
of  such  artifacts  is  perhaps  the  absence  of  any  pre-crystallization  upswings 
in  the  case  of  emulsions  of  binary  solutions  in  which  the  anomalous  character 
of  water  has  been  "quenched"  by  the  second  component.  The  behavior  of  such 

solutions  is  systematically  presented  for  the  case  of  +  H20,  in  a  pre- 

(9) 

vious  paper. 

The  most  serious  source  of  error  in  these  measurements  at  low  temperature 
would  appear  to  be  in  the  temperature  calibration.  This  must  be  assumed  the 
same  for  downscans  as  it  is  for  upscans  because  due,  to  supercooling  of  standards, 
calibrations  are  always  conducted  on  upscans. 


RESULTS 

The  new  data  for  water  are  collected  in  Table  1,  and  displayed  graphically  in 
Fig.  1.  Discussion  of  the  form  of  this  heat  capacity  function,  and  the  manner  in 
which  it  is  affected  by  changes  of  external  pressure,  and  additions  of  second 
components,  is  given  elsewhere.^ 

A  redetermination  of  for  D20  has  reproduced  earlier  reported  DSC  data 
to  better  than  the  experimental  uncertainty  of  ±2%.  The  new  data  for  bulk  and 
emulsion  samples  are  compared  with  the  old,  and  with  literature  data  for  bulk 
samples,  in  Table  2. 


(1) 


DISCUSSION 

Anomalous  increases  in  liquid  heat  capacity  at  low  temperature  of  the  form 
of  Fig.  1  are  not  found  in  other  molecular  liquids,  even  those  with  otherwise 
water-like  properties,  such  as  H202  and  N^,  nor  in  mixtures  of  water  with  more 
than  20%  of  these  liquids  (even  though  some  of  these  mixtures  can  be  studied 
during  continuous  cooling  down  to  -140®C  where  they  vitrify). ^ 


*~36ifc,v— '■ 
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On  the  other  hand,  all  molecular  liquids  which  are  internally  stable  against 
decomposition  must  exhibit  anomalies  of  this  character  on  sufficient  super¬ 
heating.  This  follows  from  the  general  form  of  the  fluid  equation  of  state 
which  shows  that  a  region  of  mechanical  instability  must  exist  in  the  isobaric 
V-T  diagram  for  all  pressures  below  the  critical  pressure.  Such  mechanical 
instabilities  are  generally  illustrated  by  reference  to  the  impossible  negative 
compressibility  region  in  the  ordinary  van  der  Waals  equation  P-V  isotherm 
for  T  <  Tc<  For  qualitative  purposes,  we  illustrate  the  corresponding  T-V 
behavior  in  Fig.  2  for  a  series  of  isobars  at  P  <  calculated  for  water  using 
a  van  der  Waals  equation  with  parameters  obtained,  as  usual,  from  the  Tc  and  Pc 
values.  The  maxima  correspond  to  minima  in  the  corresponding  P-V  isotherms, 
hence  to  points  at  which  the  superheated  liquid  becomes  mechanically  unstable, 
and  must  spontaneously  rupture.  Of  course,  the  simple  van  der  Waals  equation 
is  inadequate  to  describe  fluid  water  in  detail,  and  it  is  therefore  not  sur- 
pring  that  T-V  maximum,  hence  superheating  limit,  of  277°C  predicted  by  this 
equation  for  1  atm  pressure,  is  less  than  the  experimental  superheating  limit, 

Tfi,  of  279. 5°C,  observed  for  small  droplets  by  different  authors5^’  or 
302°C  obtained  by  Skripov  and  coworkers  from  explosive  boiling  nucleation  rate 
studies.  The  latter  value  is  close  to  that  calculated  by  Apfel  for  nuclea- 

tion  in  bulk.  A  The  theoretical  limit  (the  spinodal  temperature  Tg) ,  should 
actually  lie  somewhat  above  the  homogeneous  nucleation  temperature  (which  is 
the  kinetically  limited  reflection  of  the  impending  catastrophe)  and  has  been 

estimated  at  305°C  by  Eberhart  and  Schnyders^15 ^  from  equation  of  state  cal 

(14) 

culations,  and  at  324°C  at  1  atm  by  Skripov.  In  this  region  the  impending 

instability  must  be  anticipated  by  large  fluctuations  in  density  and  entropy 
and  these  will  be  reflected  in  the  magnitudes  of  the  related  thermodynamic 
properties,  compressibility  and  heat  capacity,  in  which  we  are  interested 


here. 
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Although  direct  measurements  of  heat  capacity  or  expansivity  in  this 
region  have  not  yet  been  made  they  may  be  calculated  using  the  equation  of 
state  described  by  Haar  et  al.^^  The  compressibility  calculations  up  to 
240°C  from  the  same  equation  are  in  very  good  agreement  with  values  obtained 
recently  by  Hareng  and  Leblond^^  by  combining  direct  light  scattering  measure¬ 
ments  of  adiabatic  compressibility  with  extrapolated  density  data.  Values  of 

/ 18) 

C  calculated  by  Haar'  at  normal  pressure  are  plotted  in  Fig.  3  together 
P 

with  the  low  temperature  data.  It  is  notable  that  even  the  upswing  in  at 
low  temperatures  is  predicted  by  this  equation.  The  increases  at  high  tempera¬ 
ture  are  more  gradual  than  those  predicted  using  a  scaling  equation  suggested 
by  Lienhar d .  ^ ^ 

The  compressibility  data  of  Hareng  and  Leblond,^^  which  were  in  good 

agreement  with  earlier  data  to  180°C  of  Trinh  and  Apfel^  and  more  recent 

(14) 

data  to  240°C  of  Skripov  and  coworkers,  were  shown  to  fit  very  well  an 
equation  of  the  familiar  critical  law  form 

<T  =  A(T/Tg  -  1)"Y  (3) 

with  T  =  588  K  and  y  =  1.  This  form  had  also  been  tested  by  one  of  the 
s 

(3) 

present  authors  on  calculated  heat  capacities  based  on  Lienhard's  equation, 

under  the  impression  that  a  value  of  y  =  1  was  predicted  by  mean  field  theory. 

However,  it  now  seems  that  y  =  1  is  only  applicable  at  the  critical  density. 

Away  from  the  critical  density,  i.e.  approaching  a  spinodal  point,  it  appears 

from  theoretical  analyses  by  Compagne^^  and  Speedy that  the  exponent 

should  be  1/2  for  ic„,  and  also  for  a  and  C  . 

T  p 


If  we  take  the  calculated  values  of  Haar  covering  the  same  Cp  range  as 
has  been  measured  at  low  temperatures,  we  find  that  after  subtracting  a  con¬ 
stant  "background"  contribution  of  75  J  K  ^  mol  ^  the  excess  Cp  can  be  fitted 
to  Eq.  (1)  with  Tg  =  600,0  K  and  y  -  1.0  -values  similar  to  those  found  by 
Hareng  and  Leblond-as  shown  in  the  insert  to  Fig.  3.  On  the  other  hand, 

extension  of  the  equation  of  state  calculations  (into  regions  where  no  validity 

( 14) 

would  be  claimed  by  its  authors)  shews  C  -*•  00  at  595.9  K,  cf  597  K  by  Skripov 

P 

and  the  limit  is  approached  with  a  smaller  Eq.  (1)  slope,  y  z  1/2,  as  the  new 

theoretical  treatments  suggest.  It  is  unfortunate  that,  because  of  metastability 

and  nucleation  kinetics,  it  will  not  be  possible  to  resolve  the  problem  of  the 

appropriate  power  law  by  experiment,  as  in  the  case  of  critical  phenomena. 

(9) 

It  has  been  shown  elsewhere  that  the  excess  heat  capacity  of  supercooled 
water  obtained  by  subtracting,  from  the  Table  1  values,  a  background  Cp  (deter¬ 
mined  from  binary  molecular  and  ionic  solution  data  described  elsewhere^ ^ )  is 
also  well  described  by  Eq.  (1)  with  a  Tg  value  of  226  K  and  y  ~  1.0.  Interest¬ 
ingly  enough,  this  singular  temperature  is  just  5  degrees  below  the  lowest 

(21 

freezing  point  ever  achieved  for  pure  water  samples  (cloud  chamber  experiments). 

The  implied  possibility  that  the  low  temperature  phenomenon  might  be  due  to  an 

(22) 

instability  was  discussed  in  some  detail  in  an  earlier  paper  and  has  now 

been  incorporated  into  broad  unified  treatment  of  high  temperature,  low  tem- 

(23) 

perature  and  negative  pressure  instabilities  in  liquid  water,  by  Speedy. 

For  the  emulsion  droplets  studied  here,  which  freeze  commencing  at  -38.5°C, 
e  (*  T/T  -1)  at  the  phase  change  is  0.038.  The  corresponding  value  for  e  at 
the  superheating  limit  of  liquid-suspended  droplets  is  0.078  if  we  use  Tg  =  600  K 
from  Haar's  equation  and  T^  =  553  from  droplet  superheating  limit  studies . ^ ^ 
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However,  this  value  of  T  may  be  too  low  since  it  corresponds  with  nucleation 

H 


at  the  silicone  oil  +  water  surface. 


(12) 


For  bulk  nucleation  T,  is  estimated 

h 


to  be  about  298°C  (571  K)  by  Apfel,  and  explosive  nucleation  measured  by  Skripov 

et  al^^  yields  =  573  K,  for  which  case  e  =  0.045,  very  close  to  the  value 

at  freezing.  On  the  other  hand  the  liquid  heat  capacities  at  the  droplet 

superheating  and  supercooling  limits  are  remarkably  similar.  This  could  be 

rationalized  by  arguing  that,  since  is  proportional  to  the  mean  square 

entropy  fluctuations,  and  since  a  critical  nucleus  is  essentially  a  special 

entropy  fluctuation,  a  correlation  of  nucleation  probability  and  C  might  be 

P 

expected. 

A  quantitative  theory  for  the  low  temperature  heat  capacity  behavior,  in 

particular  the  apparent  accord  with  a  critical  law,  is  currently  lacking. 

(21) 

Stanley  and  Teixeira  have  recently  shown  that  the  singularity,  if  real, 
could  be  associated  with  a  percolation  threshold  for  clusters  of  "four-bonded" 
water  molecules,  hence  amorphous  ice-like  density  fluctuations,  in  the  liquid. 

On  the  other  hand.  Speedy has  now  shown  that  all  the  volumetric  observations 
on  superheated,  stretched  (negative  pressure),  normal  and  supercooled  water  can  be 
self-consistently  explained  by  postulating  the  existence  of  a  spinodal  (k^,  -*  <») 
line  bounding  the  mechanically  stable  liquid  regime,  and  deriving,  by  a  Landau 
expansion,  an  equation  of  state  which  yields  Pg  on  fitting  to  the  available 
volumetric  data.  The  form  of  this  Pg,  Tg  boundary  line  is  that  of  a  normal 
liquid  at  high  temperature,  but  at  lower  temperature  and  negative  pressures, 
it  departs  from  the  normal  van  der  Waals  behavior  (which  is  to  constantly 
increasing  negative  Pg  with  decreasing  Tg)  and  at  Tg  =  35°C  Pg  *  -2 IDO  bar, 
it  reverses  slope.  The  reversal  point  corresponds  with  the  intersection  of 
the  spinodal  with  the  extrapolation,  to  negative  pressures,  of  the  observed 
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temperatures  of  density  maxima.  The  spinodal  (instability)  thus  reappears  at 
positive  pressures  at  temperature  below  the  freezing  point.  The  peculiar 
reversal  would  be  a  feature  common  to  liquids  in  which  the  lowest  energy  state, 
which  is  approached  with  decreasing  temperature,  is  correlated  with  the  largest 
volume  state.  Unfortunately,  there  are  very  few  such  liquids  (SiO^  and  some 
tetrahedrally  coordinated  antimonides  and  tellurides)  so  the  phenomenon  cannot 
be  broadly  studied. 

Both  the  percolation  and  the  spinodal  models  have  attractive  and  challenging 

features,  and  will  doubtless  stimulate  increasing  experimental  and  theoretical 

work  in  this  area.  However,  a  tendency  evident  in  some  of  the  current  literature 

to  see  these  two  approaches  as  somehow  mutually  exclusive  should  be  discouraged. 

The  language  used  in  ref.  22  to  describe  the  molecular  nature  of  the  increasing 

density  fluctuations  which  the  power  law  for  the  compressibility  required  to 

be  present,  ("density  fluctuations  presumably  arise  and  decay  by  the  migration 

of  large  0***H-0  angle  (-*  180°)  hydrogen  bonds  which  'prop  open'  the  expanded 

structure  when  present  in  high  local  concentration  but  allow  collapsed  dense 

regions  to  form  when  local  concentrations  are  smaller",)  accurately  previews 

(24) 

the  topological  picture  provided  by  the  percolation  model.  Future 

refinements  of  the  current  simplified  percolation  model  to  incorporate  the 
actual  continuous  distribution  of  hydrogen  bond  energies  and  other  interaction 
terms  may  well  lead  to  an  equation  of  state  containing  an  instability  at  low 
temperatures,  as  indeed  i£  contained  in  the  optimized  empirical  equation  of 

state  of  Haar  and  colleagues  (which  we  have  only  used  in  this  paper  to  discuss 

(25) 

the  high  temperature  anomalies  ).  In  this  case  the  two  approaches  would 
obviously  be  mutually  reinforcing. 

Both  approaches  are,  however,  antagonistic  to  heterophase  fluctuation 
interpretations  of  water  anomalies  since  the  latter  contain  no  reason  for 
excluding  any  liquids  susceptible  to  crystallization  from  showing  qualitatively 


similar  anomalies. 
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Heat  Capacity  of  Water  in  J  K 


(a)  Handbook  of  Physics  and  Chemistry,  CRC  press,  58th  Edition. 


Temperature 


T/K  T/°C 


298.15 

295 

290 

289 

285 

284 

280 

279 

275 

270 

265 

263 

260 

255 

253 

250 

245 

242 

240 


Table  2 


Heat  Capacity  of  D2O  in  J  K  ^  mol  1 


Investigator  I 


lit. 


Run  I  Run  II 


83.75 

84.1 


84. 3+ 
82.75+ 
82.4* 
83. 58+ 

83.62+ 

84 . 36+ 


emulsion 

Run  I 

Run  II 

85.4 

83.5 

84.9 

83.7 

85.1 

84.1 

86.0 

85.0 

87.2 

86.5 

89.6 

88.7 

93.3 

92.3 

99.9 

99.2 

109.3 

109.5 

119.6 

129.3 

ref .  1 


emulsion 


92.1 

95.5 

101.3 

110.0 

129.2 


~  Termicheskie  Konstanty  Westchestv  (Thermal  Constants  of  Substances)  Academy 
of  Sciences,  USSR  Moscow,  1965,  Part  I,  p.  33. 

*  NBS  Circular  500 

+  Kemp  and  Roth,  J.  Am.  Chem.  Soc.  1936,  ^8,  1829. 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Heat  capacity  of  water  in  normal  and  supercooled  temperature 
ranges  at  1  atm  pressure.  •  literature  data  for  bulk  sample 
O,  □,  A  .  data  obtained  in  this  study  by  investigator  I  for 
bulk  sample  (O)  and  emulsion  samples  using  MCR  carrier  (□) 
and  n-octane  carrier  (A):  H  data  obtained  in  this  study  by 
investigator  II  for  emulsion  sample  in  n-octane  carrier. 

Isobars  for  water  at  subcritical  pressures  according  to  the 
Van  der  Waals  equation  showing  the  metastability  limits  at 
various  pressures. 

Combination  of  low  temperature  range  data  of  this  study  with 
metastable  superheated  range  data  calculated  from  the  equation 
of  state  of  Haar  et  al  for  superheated  range.  The  metastable 
state  low  and  high  temperature  limits  observed  for  small  droplets 
are  -42°C  and  280°C,  respectively  (marked  T^) . 

Insert:  Test  of  critical  type  law  for  high  temperature  equation- 
of-state  0^  values  after  subtracting  a  constant  "background"  of 
75  cal  mol  ^  deg  The  plot  uses  Tg  =  600.0  K,  and  yields  a 
slope  of  unity  as  indicated  by  the  solid  line  in  the  inset. 

Use  of  (calculated)  C  values  beyond  the  C  range  of  those  measured 

p  — i -  p 

in  low  temperature  experiments,  yields  slightly  lower  Tg 
(Tg  -*•  596  K)  and  smaller  slopes  (y  -*■  0.4). 


TYPE  OF  DATA  INVESTIGATOR 

*  * '  BULK  SAMPLE 

°  °  °  MCH  CARRIER  I 

*  *  “  N-OCTANE  CARRIER. 

*  •  ■  N-OCTANE  CARRIER  0 

*  *  *  LITERATURE  DATA 

- PREVIOUS  AVERAGE  (DSC- 1  and 

DSC-2) 

- PREVIOUS  DSC-2 


I  1% 

V, 

n  rf  *  ?  . 
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